Effect of the Centrifugal Force on Domain Chaos in Rayleigh-Benard convection 
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Experiments and simulations from a variety of sample sizes indicated that the centrifugal force 
significantly affects rotating Rayleigh-Benard convection-patterns. In a large-aspect-ratio sample, 
we observed a hybrid state consisting of domain chaos close to the sample center, surrounded by an 
annulus of nearly-stationary nearly-radial rolls populated by occasional defects reminiscent of undu- 
lation chaos. Although the Coriolis force is responsible for domain chaos, by comparing experiment 
and simulation we show that the centrifugal force is responsible for the radial rolls. Furthermore, 
simulations of the Boussinesq equations for smaller aspect ratios neglecting the centrifugal force 
yielded a domain precession-frequency f ^ e'^ with ~ 1 as predicted by the amplitude-equation 
model for domain chaos, but contradicted by previous experiment. Additionally the simulations 
gave a domain size that was larger than in the experiment. When the centrifugal force was included 
in the simulation, /i and the domain size closely agreed with experiment. 

PACS numbers: 47.54.-|-r,47.32.-y,47.52.-|-j 



I. INTRODUCTION 

Rayleigh-Benard convection (RBC) occurs when the 
temperature difference AT across a horizontal layer of 
fluid confined between two parallel plates and heated 
from below reaches a critical value ATc P, ■ Rotat- 
ing the sample about a vertical axis at a rate $7 > fic 
induces a state of spatio-temporal chaos 0, 1^ known as 
domain chaos. 

Domain chaos featurespatches of straight rolls oriented 
at various angles HHlElSlilllSllflllllll- The 
orientation of the rolls fluctuates in time and space ap- 
proximately at discrete angles due to the Coriolis force 
which induces the Kiippers-Lortz instability. The size 




Figure 1: Images of domain chaos in Rayleigh-Benard con- 
vection with r = 36, e = AT /ATc - 1 = 0.05, and Q = 16.25. 
Left: the temperature profile at the mid-plane from simula- 
tion of the Boussinesq equations with no centrifugal force. 
Center: shadowgraph image from experiment. A movie of 
this sample at 150 x actual speed is available [l^ . Right: the 
temperature profile at the mid-plane from simulation of the 
Boussinesq equations including centrifugal force. 
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and shape of the domains also fluctuates and moving de- 
fects pepper the entire pattern creating a state of persis- 
tent chaos in both space and time. Figurenshows several 
snapshots of domain chaos. 

The theoretical models that describe domain chaos 
usually neglect the centrifugal force, relying on the Corio- 
lis force as the dominant effect brought by rotation. Even 
when the full Boussinesq equations of motion are con- 
sidered, the centrifugal force often is evaluated on the 
assumption that the density is constant throughout the 
sample Q. In that case it has no influence on the neu- 
tral curve and on the pattern that forms above it. In 
this paper we show, both from experiment and from nu- 
merical simulations of the Boussinesq equations, that the 
centrifugal force does have a significant influence on the 
quantitative features of domain chaos. This is so even 
for modest aspect ratios T = ro/d, where vq and d are 
the radius and height of the convection sample. When T 
becomes large enough the centrifugal force qualitatively 
alters the observed patterns. In order to reproduce these 
features in the simulation, it is necessary to consider the 
temperature dependence of the density, as is done explic- 
itly in the Appendix. 

Qualitatively the influence of the centrifugal force is 
readily understood. The vertical density gradient, result- 
ing from the imposed temperature difference, induces a 
radial large-scale circulation (LSC). The wave director 
of the RBC rolls tends to align orthogonally to the LSC 
[T^ . This tendency competes with the Kiippers-Lortz 
mechanism that tends to create disordered and fluctuat- 
ing domains. 

For samples with F up to about 40 we observed pat- 
terns in the experiment that were regarded as typical 
of domain chaos and similar to those fou nd p reviously 
^3,0,0- However, as observed before |l2j |. the pat- 
terns contained domains that were signiflcantly smaller 
than those of simulations of the Boussinesq equations. 
At larger F = 80 we found experimentally a hybrid state 
where the centrifugal force was strong enough relative to 
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the Coriolis force to qualitatively affect the pattern as 
shown in Fig. |21 There the central section looked like do- 
main chaos, but the annular region along the perimeter 
was mostly made up of nearly-stationary nearly-radial 
rolls. Defects glided azimuthally across the radial rolls 
(a movie of this state can be found in the online ver- 
sion) . Although simulations of the Boussinesq equations 
for such large T were not feasible, the inclusion of an en- 
hanced centrifugal force in simulations for F = 20 showed 
a similar radial roll structure near the sample boundary. 

The outer region of the hybrid state had much in com- 
mon with undulation chaos [isl IT^ IT^ observed in in- 
clined Rayleigh-Benard convection. In the inclined sys- 
tem there is a component of gravity in the plane of the 
sample. The undulation chaos consists of defects gliding 
across straight rolls, with the roll axes aligned parallel to 
this in-plane component. Similarly, there is an alignment 
of the roll axes in the rotating case that is nearly paral- 
lel to the (radial) centrifugal force, with defects gliding 
azimuthally across the rolls. The resemblance between 
these phenomena exists in spite of the fact that these 
forces possess a very different character because the cen- 
trifugal force depends on radial position while gravity is 
uniform throughout the sample. 

In the remainder of this paper we first give some ex- 
perimental and theoretical/numerical details, and then 
present several results of quantitative pattern analysis 
that characterize the extent of the influence of the cen- 
trifugal force for 15 < fi < 19 (fi = lo(P /v where lo is the 
angular rotation rate and v the kinematic viscosity) and 
e = AT/ATc — 1 0.5. On the basis of this analysis we 
show for this parameter range that the crossover from 
the Coriolis-force-dominated central region of domain 
chaos to centrifugal-force-dominated near-radial rolls oc- 
curs near r/d = 35 (r is the radial coordinate) when 
e ^ 0.1, and at somewhat smaller r/d for smaller e. We 
conclude that previous samples had been too small to 
observe the qualitative features of the hybrid state. Un- 
fortunately we also have to conclude that the experimen- 
tal study of the Coriolis-force effect in compressed gases, 
where large aspect ratios are attainable, will be severely 
contaminated by centrifugal-force effects when F > 40. 



II. DETAILS OF EXPERIMENT AND 
SIMULATION 

The sample exhibiting the large-F hybrid state (see 
Fig. El had a height d = 553 ^m with F = 80. It 
contained sulfur hexa-fluoride (SFg) at a pressure of 
25.44 bars and a mean temperature of 43.00°C where 
the Prandtl number a = v/k was 0.88 {k, is the ther- 
mal diffusivity) . We used the shadowgraph technique to 
image the convection patterns in an apparatus described 
elsewhere [Tsl. Although Fig. H shows an image taken at 
f2 = 16.25, we observed a similar phenomenon over the 
entire range 15 < VI < 19 of accessible rotation rates that 
were substantially above flc- 




Figure 2: (Color online) Hybrid state of domain chaos with 
radial rolls and defects along the outer annulus. System pa- 
rameters are F = 80, e — 0.05, and fl — 16.25. The grey-scale 
(false color in the electronic version) overlay is 6 (r), the an- 
gular component of the local wave-director field. A movie is 
available [lyl showing the dynamics of this pattern at 150 x 
actual speed. 



We also studied a medium and a small sample exper- 
imentally that did not exhibit the hybrid state. The 
medium sample had F — 61.5, a thickness of 720/im, 
a pressure of 20.00 bars, a mean temperature of 38.00°C, 
and a = 0.87. The small sample had F = 36, was 1230/im 
thick, operating at 12.34 bars, with a mean temperature 
of 38.00°C, and a = 0.82. 

In addition to the experiment, we also simulated the 
Boussinesq equations to study the system theoretically. 
A discussion of the relevant equations of motion is given 
in the Appendix. Details of the numerical code are de- 
scribed elsewhere ^2d,^2l]. No-slip velocity boundary con- 
ditions and conducting lateral thermal boundary condi- 
tions were utilized. We used a spatial resolution of 0.1 
and a time resolution of 0.005. In most simulations that 
included the centrifugal force, the centrifugal term was 
made larger than would be physically realistic, in the 
present experimental fluid, in order to model the effect 
of a larger aspect ratio. 

Two different configurations were used to compute the 
precession frequency /. In one case, which included 
both the centrifugal and Coriolis force, we used F = 20, 
a — 0.93, n — 17.6, and the centrifugal term was twice 
the physically realistic value in order to model a F that 
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would be twice as large. In the other case, where only 
the Coriohs force was included we used F = 40, 
a = 0.93, and = 17.6. We compare the results of these 
simulations to the experimental results from Ref. |lOj| . 
which also used T = 40, cr = 0.93, and f2 = 17.6. 

We also ran some simulations with the exact parame- 
ters of the r = 36 experimental sample for direct com- 
parison. Finally we ran several additional cases for the 
r = 20 simulation using a centrifugal term with various 
strengths. 



III. THEORETICAL CONSIDERATIONS 



Althoui 
iii 



revious theoretical work on domain chaos 
III22, 23, 24, 25] neglected the centrifugal 



force because the Froude number Fr 



r/g (g is the 



acceleration of gravity) is small in typical experiments, 
other research |2y, |23| indicated that the Froude number 
may not be the only relevant parameter. Hart pro- 
vided a numerical estimation of the parameter regimes 
where the centrifugal force is relevant. He made the ap- 
proximation that n is large (> 500), so we cannot directly 
use that estimate because the present values are only 
about 20. We instead consider 



A 



2uoei/2 



(1) 



(/3 = aATc where a is the isobaric thermal expansion 
coefficient). A derivation of the relevant centrifugal and 
Coriolis terms is given in the Appendix. This parame- 
ter A is the small-£ approximation to the ratio of the 
magnitude of the maximum of the centrifugal term, eval- 
uated near the outer edge of the sample, to the magni- 
tude of the Coriolis term evaluated where it reaches a 
maximum in the horizontal direction, z — zq- (Note that 
the temperature of the conduction profile is equal to — z 
for our system). The quantity A is an indicator of the 
transition from domain chaos to a hybrid state. We can 
obtain an approximate numerical value for uq by per- 
forming a linear stability analysis to get the functional 
form of the velocity and by using the simulation to get 
the normalization. For the case of F = 40, a = 0.93 
and f2 — 17.6, we find uq = 12.3 at zq = 0.3. These 
values should also approximately apply for our experi- 
mental parameters, which do not vary significantly from 
the parameters used to obtain the numerical values, ex- 
cept for r whose dependence is explicit. For the F = 80 
sample, A — 0.61 for e ~ 0.05, indicating that the cen- 
trifugal force is almost as influential as the Coriolis force, 
but for the F = 36 sample and e = 0.05, A — 0.12 indi- 
cating that the Coriolis force dominates. Both samples 
have similar values of and a, so it is not surprising that 
we observe a F-dependent transition to a hybrid state in- 
duced by competition between the centrifugal force and 
the Coriolis force. 




Figure 3: Simulation of the Boussinesq equations with F — 
20, (T — 0.93, Q = 17.6, e = 0.055, and an artificially large 
centrifugal force in order to model the efltect of an inaccessibly 
large F. Left: zero centrifugal force. Center: centrifugal force 
is 4 times the physical value. Right: centrifugal force is 10 
times the physical value. A movie of this sample at 100 x 
actual speed is available |2^. 



IV. RESULTS 

We qualitatively tested our hypothesis that the cen- 
trifugal force is responsible for the hybrid state by sim- 
ulating the Boussinesq equations with centrifugal force 
included. It was not possible to reach F = 80 in the 
simulation in order to directly reproduce Fig. [3 Instead 
we ran the simulation at F = 20 and with the centrifu- 
gal force artificially large in order to model the stronger 
centrifugal force at larger F. Some results are shown in 
Fig. |3| The qualitative agreement between Fig. |2 and 
the right image of Fig. O is striking. Both exhibit do- 
main chaos at the center but a radial roll structure in an 
annulus near the side wall. 

In order to make a direct comparison with the exper- 
iment we also simulated the Boussinesq equations for 
F = 36 at e = 0.05 with the exact parameters of the 
F = 36 experimental sample, both with and without the 
centrifugal force included. FigureHcompares examples of 
the resulting pattern between simulation, with and with- 
out the effect of centrifugal force, and experiment. As 
noted before |0, we observed that the domain size was 
too large when we neglected the centrifugal force as seen 
by comparing Figs.^ and^. The domain size in pat- 
terns where we included the centrifugal force, as shown 
in Fig. qualitatively matched the experimentally ob- 
served shadowgraph image in Fig. ^p. 

Recent experiments ^3 disagreed with the theoretical 
prediction H, ^ ^ e~^^'^ for the scaling of the correla- 
tion length ^ and f ^ with /J. = 1 for the scaling of the 
domain precession frequency /. We do not address the 
^ scaling in the present work because finite-size effects 
cloud the issue .22. .29<l. However, fj, is unaffected by the 
finite-size effects~ |22L |2 9| . The method of determining / 
is discussed in Refs. pXll2l| . We find that / from the ex- 
periment agrees with the simulation when the centrifugal 
force is included. The scaling for / is shown in Fig. 21 for 
simulations both with and without the centrifugal force. 
The simulations that neglected the centrifugal force gave 
a scaling exponent /i = 1.15, in good agreement with the 
prediction of /i = 1 from amplitude equations 8'] which 
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Figure 4: Effect of the centrifugal force on the scaling of the 
precession frequency / with e. All data are from samples with 
a = 0.93 and Q — 17.6. Circles: from simulations that include 
both the centrifugal force and the Coriolis force, with F = 20 
and a centrifugal term that is twice as large as it would be 
in the experimental fluid. Diamonds: from simulations that 
include only the Coriolis force, with F = 40. Squares: data 
from previous experiments, with F = 40 fToj . Solid lines: 
least-squares fits of / — Ce'^ to the data. 



also neglected the centrifugal force. There is also good 
agreement between the experimental data and the results 
of the simulation where both the centrifugal and Coriolis 
force were included. Not only did the numerical simula- 
tion yield a scaling exponent of /i = 0.62, very close to 
/i = 0.58 from the experiment, but it also reproduced the 
actual values of / in the experiment remarkably well. 

Aside from the dramatic effect on the pattern seen in 
Fig. 121 in the experiment we also observed a downward 
shift by three or four percent of the critical Rayleigh num- 
ber Rr. = aq /!^Tr.(P / Kv for the large sample. Previous 
work |26l |27| showed that the centrifugal force may ei- 
ther stabilize or destabilize the convective flow depending 
on the parameter regime, although the range of param- 
eters specifically corresponding to the present case was 
not directly investigated. However it was suggested that 
for small the centrifugal force lowers Rc [23 ■ 

We did not observe any radial dependence of Rc for 
the r = 80 sample. We measured Rc in both a central 
circular region of half the sample radius and in the outer 
annulus surrounding that region. We accomplished this 
by computing the time averaged total power in the spec- 
ified region of the experimental images as a function of 
AT. The total power is proportional to A/" — 1 where 
TV is the Nusselt number. Through a quadratic fit we 
extrapolated to the background value of the total power 
in order to measure ATc, from which we calculated Rc 
from knowledge of the fluid properties. The values for 
Rc in both regions agreed within the statistical error of 
the data, with the biggest difference being a bit less than 
0.2%, i.e. an order of magnitude less that the down shift 
of Rc. 
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Figure 5: Shifted critical Rayleigh number due to the centrifu- 
gal force. Solid line: prediction from linear stability analysis. 
Circles: measured from the F = 61.5 sample. Squares: mea- 
sured from the F = 80 sample. 

The linear stability analysis P| for this system was 
carried out for a laterally infinite system using the Cori- 
olis force only and neglecting the centrifugal force. For 
r — 61.5, our measured Rc agreed well with this analysis 
as is seen in Fig. |S1 Similarly good agreement had been 
found earlier for F = 40 This is interesting because 
we do observe other centrifugal effects for F = 61.5, as 
discussed below. For F = 80 there was a clear downward 
shift of Rc, relative to the theory and to experiments for 
smaller F, of roughly 3%. This trend with F indicates 
that the disagreement with theory cannot be attributed 
to the finite size of the experimental systems because any 
finite-size effect on the onset should have been larger at 
smaller F. Thus it suggests that for F = 80 the centrifu- 
gal force, which increases with increasing F, has become 
strong enough to affect the onset. 

We quantified our analysis of the experimental im- 
ages by utilizing the angular component (r) of the lo- 
cal wave-director field computed with an algorithm de- 
scribed in Ref. but using higher angular resolution. 
Figure|21shows an example of (r) overlayed on top of the 
shadowgraph pattern in false color. We used the variance 
field of 6 (r), from a time series of consecutive images, to 
address the issue of which region of the pattern exhibits 
chaotic dynamics and which region is mostly dominated 
by near-stationary radial rolls. Chaotic patterns exhibit 
a large variance because the roll orientation fluctuates 
constantly, but stationary rolls yield a small variance. 

The mean angular sum fleld is 

N 

5(r) = ^e2^''"W (2) 

n 

where n is the image index and A'^ is the number of im- 
ages. The complex exponential introduces the necessary 
periodicity for summing angles and the factor of 2 treats 
the field as a director field instead of as a vector field. 
The mean angular field is 
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Figure 6: The variance field (O"^^ (r) for several system sizes 
from experiment with e = 0.05 and Q. = 16.25. Black corre- 
sponds to {O^^ (r) = 2 and white corresponds to {d'^^ (r) = 0. 
Left: r = 36. Middle: V = 61.5. Right: T = 80 




, 1 /lmS'\ 



(3) 



The angular variance field is 



AT 

(02^ (r) ^ ^^|e2.e„(r) _g2.fl(r) 



1 



(4) 



Simple algebra shows that < (0^) (r) < 2 where a value 
of indicates that the domain orientation is stationary 
while a value of 2 indicates that it is maximally fluc- 
tuating. The variance field lO"^) (r) is approximately az- 
imuthally symmetric for several investigated system sizes 
as shown in Fig.|^ There one sees for the larger two sam- 
ples that the variance is large near the center and smaller 
near the side wall. 

FigureHdisplays (6*^) (r), the azimuthal average of the 
angular variance field, for several system sizes. For the 
smallest experimental sample and for a simulation with 
matching parameters, both with F = 36, the large vari- 
ance across the entire sample revealed the presence of 
domain chaos throughout. However, as seen in Fig. ^ 
even in this case the centrifugal force has reduced the 
domain size. Likewise, the central region of the largest 
sample, with F — 80, exhibited a large variance. Domain 
chaos dominated near this central region, but along the 
perimeter of the sample the variance was very small in- 
dicating coexistence with a more nearly stationary pat- 
tern. The variance for both the F = 61.5 and F = 80 
samples showed a clear crossover from domain chaos to 
radial rolls, however the F = 61.5 sample did not fully 
transition out of domain chaos even at the edge of the 
sample. 

We also measured the size of the domain chaos region 
of the hybrid state by finding the width h of lO'^) (r) 
chosen such that (9^) (6) = (6^) (0) /2. Figure El shows 
the £ dependence of b. The behavior of A (see Eq. at 
small s agrees with such e dependence because of the fac- 
tor e^/^ in the denominator, associated with the Coriolis 
force, which gives A ^ e'^l"^ for small e. The data in 
Fig. El indicate that, for the present sample parameters, 



Figure 7: Azimuthal average {O^^ (r) of the angular variance 
field (shown in Fig. |SJ for various system sizes from both ex- 
periment and simulation. All data were taken with e = 0.05 
and Q, = 16.25. Solid line: experiment, F = 36. Dashed line: 
experiment, F — 61.5. Dotted line: experiment, F = 80. 
Dashed-dotted line: simulation including centrifugal force, 
F = 36. 




Figure 8: Width b of (O^) (r) for F = 80 at various e. Crosses: 
!^ = 15. Circles: = 16.25. Squares: = 17.7. Pluses: 
Q, = 19. 



the centrifugal force dominates at small e while the Cori- 
olis force becomes more important at large e, although 
for large enough F neither are negligible. 

In addition to looking at the fluctuations in the wave- 
director angle, we also measured the angle of the wave 
director relative to the angle 9f of the side-wall normal 
by computing a time-averaged dot-product-like quantity 



1 ^ 



(5) 



Here we are using Of as a reference for the orientation of 
the centrifugal force. The quantity D (r) is 1 when the 
wave director is parallel to the sidewall normal and -1 
when it is perpendicular. 

In the presence of the LSC induced by the centrifugal 
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Figure 9: The time-averaged dot product D (r) (the azimutlial 
average of D {r)) for the same data as shown in Fig. All 
data were taken with e — 0.05 and = 16.25. Solid line: 
experiment, F = 36. Dashed line: experiment, F = 61.5. 
Dotted line: experiment, F — 80. Dashed-dotted line: simu- 
lation including centrifugal force, F = 36. 



-0.25 





Figure 10: Average minimum D (r) for F = 80 at various e. 
Crosses: SI — 15. Circles: SI = 16.25. Squares: — 17.7. 
Pluses: Q = 19. 



Figure 11: Azimuthal average of the auto-correlation function 
of the wave-director angle field for F = 36, e — 0.05, Q — 
16.25, a — 0.82. Solid line: from simulation of Boussinesq 
equations with centrifugal force included. Dashed line: from 
simulation but with centrifugal force neglected. Dotted line: 
from experiment. 



statistical fluctuation in the D (r) curves, the plateau in 
those curves was averaged over a small radial length of 
roughly bd to yield an average minimum value for D (r) , 
min[Z?(r)], at the given e and fl. As a function of e, 
min [Z? (r)] exhibited a minimum, between e = 0.1 and 
£ = 0.3 depending on f2, which corresponded to the near- 
est realization of orthogonality between the sidewall nor- 
mal and the wave-director. For small e the rolls started 
to turn in and become nearer to concentric rings along 
the outer annulus. Such behavior is indicative of the e 
dependence of the relative strengths of the centrifugal 
and Coriolis forces. 

We also compared auto-correlation functions of the 
wave-director field between the experiment and simula- 
tion. Using 9 (r) we computed 



force, we expect the wave director to line up nearly per- 
pendicular to that flow [l^ . Some deviation from this 
might be expected due to the action of the Coriolis force 
on the Rayleigh-Benard rolls. The LSC itself would be 
expected to nearly align with the centrifugal force and 
thus should be nearly orthogonal to the side wall. How- 
ever, here also some deviation from orthogonality to the 
walls would be expected due to the Coriolis-force influ- 
ence. As expected, Fig.l^shows that, as F increases, the 
wave director along the edge of the sample indeed does 
become more nearly, but not perfectly, perpendicular to 
the sidewall normal. For the largest F, D (r) reached a 
plateau around -0.8, which corresponded to an angle of 
about 74°, somewhat shy of a perpendicular angle. 

The angle also depended on e as indicated by Fig. 1101 
which shows the minimum value of D (r) for many curves 
similar to those shown in Fig. 1^1 To reduce the effect of 



which is slightly different from the formula used in 
Ref. . The difference is the inclusion of a factor M (Sr) , 
the number of data points integrated over to find the cor- 
relation at the coordinate 6r, which corrects for the di- 
minished quantity of data available at increasingly sep- 
arated correlations. This correction results in what is 
known as an unbiased correlation function, as opposed to 
the biased correlation fimction given when Af (Sr) = 1. 
Both approaches have merits, in particular the biased 
form must be used for computing power spectra, but we 
prefer the unbiased form in the present case. 

The angular auto-correlation function magnifies the 
difference between centrifugal force and lack of centrifu- 
gal force as shown in Fig. 1111 The agreement between the 
experiment and the full simulation with centrifugal force 
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was excellent at small and moderate Sr, and indicated a 
quantitative agreement of the domain size between exper- 
iment and simulation with centrifugal force, as observed 
qualitatively in Fig. ^ The slight deviation at large sep- 
arations should probably be disregarded due to lack of a 
sufficient quantity of data, even with the help of the un- 
biased form of the auto-correlation function to mitigate 
this issue. The simulation which neglected centrifugal 
force exhibited an overall larger correlation than the ex- 
perimental result, in agreement with the larger domains 
shown in Fig. ^ 



V. SUMMARY AND CONCLUSION 

We have shown that the centrifugal force plays an 
important role in the domain chaos state of Rayleigh- 
Benard convection. Even for relatively small F where 
A remains small throughout the sample, the centrifugal 
force affects both the pattern domain size and also the 
precession frequency. At large enough F, a hybrid state 
forms that contains both domain chaos and radial pat- 
terns reminiscent of undulation chaos. We have charac- 
terized this hybrid state using local wave-director analy- 
sis. We also observed a slight shift in Rc for this state. 

Comparison between simulations and experiment at 
moderate F indicated that the centrifugal force is respon- 
sible for the disagreement of previous experiment with 
the scaling / ~ e predicted by the amplitude model for 
domain chaos. A simulation which neglected centrifu- 
gal force agreed with that model, but a simulation which 
included centrifugal force agreed with the experiment. 
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Appendix A: ROTATIONAL CORRECTIONS TO 
THE BOUSSINESQ EQUATIONS 

We begin with the Navier-Stokes equation in a rotating 
frame and in the presence of gravity, the heat equation 
and mass conservation: 

VP 

(9( + u-V)u = h J^V^u - - 2tjz X u + cj^r, 

P 

{dt + u-V)T = kW^T, 
V-u = 0, 

(Al) 

where u is the velocity field, T is the temperature field, 
P is the pressure field. 

We apply the Boussinesq approximation, in which all 
fluid parameters are assumed to be constant except for 
the density in the buoyancy term. Unlike the standard 
application of this approximation P|, we include buoy- 
ancy from the centrifugal force as well as gravity. To 
lowest order, the density variation in these terms is 

p^p[l~a{T-f)], (A2) 

where T is the mean temperature and p is the density at 
that temperature. 

As long as a(T — T) is small, we can expand the pres- 
sure term: 

VP VP 

- — ^- — [l + a{T-T)]. (A3) 
P P 

Since the pressure is determined from a gradient, we can 
absorb the hydrostatic pressure due to the gravitational 
and the centrifugal forces into the dynamic pressure, by 
redefining P = P + pgzz — |(jj^(a;^x -I- y'^y). The terms 
proportional to the temperature will not be absorbed. 
Equation (|A1|) becomes: 

VP 

((9t + u-V)u = — + iyV^u + ga{T-T)z 

P 

- 2uj z X u- uj^a{T -f)r, (A4) 

where we have assumed the term VPa{T — T)/p\s small. 

The variables are then non-dimensionalized by speci- 
fying the length in terms of the cell height d, the tem- 
perature in terms of AT, and the time in units of the 
vertical thermal diffusion time = / k. We also define 
the Prandtl number a = v / and the Rayleigh number 
R = ag/^T(P / Kv . In addition, we define (3 = aATc, 
where ATc (and corresponding critical Rayleigh number 
Rc) is the temperature difference at which conduction 
gives way to convection. We obtain: 

a-^ {dt + \x-V)\i = -VP + V^u + R{T -T)z 

- 2nzxu-Pafl^ — (P-P)r, 
Rc 

{dt + u-V)T = V^T, 
V-u ^ 0. 

(A5) 
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The ratio of the magnitude of the centrifugal term to the 
magnitude of the Coriohs term is evaluated at the posi- 
tion where the Coriolis force reaches its maximum value, 
which is where the magnitude of horizontal velocity u±_ 
is a maximum. To first order in e^/^ we can define this 
magnitude — uo£^^^, which occurs at z = zq. At 
this point, the absolute value of T — T = zq, the temper- 



ature due to the conduction profile. (Note we have ne- 
glected convective corrections to the conduction profile, 
since these are of order e^^'^). We also want to evaluate 
the centrifugal force at the horizontal location where it 
reaches its maximum value, which is at r = F, which 
then yields Eq. Q 
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